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ABSTRACT 


Baroclinic  eddies  with  diameters  of  about  10  km  and  maximum 
current  speeds  of  about  50  cm/s  have  been  widely  observed  in  the 
central  Arctic  Ocean  north  of  Alaska  and  Canada.  The  possible 
origin  of  these  eddies  through  an  instability  of  the  mean  baro¬ 
clinic  flow  is  investigated  using  an  ocean  model  with  exponential 
profiles  of  mean  shear  and  VSisSlS  frequency.  The  model  includes 
Ekman  pumping  at  a  rigid  bottom  and  at  either  a  free  or  rigid 
upper  surface.  The  central  Arctic  Ocean  where  the  eddies  were 
found  is  baroclinically  stable  with  no  possibility  of  eddy  pro¬ 
duction.  If  the  eddies  are  spawned  by  this  mechanism,  they  must 
be  formed  at  a  site  far  from  where  they  were  observed.  On  the 
periphery  of  the  Arctic  Ocean  north  of  Alaska  the  combination  of 
greater  current  shear,  shallower  depth  and  lack  of  ice  cover 
leads  to  unstable  conditions  and  the  eddies  apparently  originate 
in  that  region.  The  instability  theory  predicts  maximum  velocity 
at  the  surface  instead  of  below  the  surface  as  observed.  Appar¬ 
ently  after  formation  in  open  water  the  eddies  are  advected 
beneath  the  ice  cover  and  dissipate  the  momentum  of  their  upper 
layers  against  the  ice.  This  is  demonstrated  by  calculations  for 
the  diffusion  of  vorticity  against  the  ice  in  the  case  of  an 
initial  exponential  profile.  A  subsurface  maximum  then  develops 
which  resembles  the  observed  profiles. 


' 
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1.  Introduction. 


The  temperature,  salinity  and  velocity  fields  of  the  oceans 
have  a  mesoscale  variability  on  scales  of  tens  to  hundreds  of  kilo¬ 
meters  horizontally  and  on  the  order  of  months  in  time.  These 
mesoscale  features  frequently  take  the  form  of  approximately  circu¬ 
lar  baroclinic  eddies  with  a  nearly  geostrophic  force  balance.  It 
is  likely  that  mesoscale  currents  play  an  important  part  in  driving 
the  large-scale  ocean  circulation  and  in  transporting  heat  as  well 
as  chemical  and  biological  parameters  on  a  global  scale.  A  review 
of  the  literature  of  mesoscale  currents  (McWilliams,  1979)  shows 
that  these  features  have  been  widely  observed  and  appear  to  be 
nearly  ubiquitous  in  the  world  ocean. 

Mesoscale  eddies  in  the  Arctic  Ocean  were  reported  by  Hunkins 
(1974)  and  Newton  et  al  (1974)  based  on  research  carried  out  from 
drifting  ice  stations.  During  the  AIDJEX  (Arctic  Ice  Dynamics 
Joint  Experiment)  project  in  1975  and  1976,  a  much  more  extensive 
survey  was  made  and  large  numbers  of  arctic  mesoscale  eddies  were 
identified  (Hunkins,  1980;  Manley,  1981).  Profiles  of  temperature, 
salinity  and  current  were  made  on  a  daily  basis  from  four  manned 
AIDJEX  ice  camps  spaced  100  km  apart  drifting  in  the  central  Arctic 
Ocean.  The  arctic  eddies  are  strongly  baroclinic  with  a  subsurface 
velocity  maximum  in  the  steep  density  gradient  between  30  and  200  m. 
Vertical  profiles  of  oceanographic  parameters  through  an  arctic 
eddy  and  in  adjacent  waters  outside  the  eddy  are  shown  in  Figs.  1 
and  2.  The  maximum  current  speed  of  58  cm/s  in  this  particular 
eddy  contrasts  with  the  geostrophic  mean  currents  in  this  region 
which  are  only  5  to  10  cm/s.  During  their  14-month  drift  which 
covered  a  total  track  length  of  8,140  km,  the  four  AIDJEX  camps  en¬ 
countered  143  eddies.  The  location  of  these  eddies  is  shown  in  Fig.  3. 
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One  of  the  eddies  was  recrossed  after  5 h  months,  demonstrating  that 
they  persist  for  at  least  that  length  of  time.  There  is  a  peak  in 
mean  vertical  kinetic  energy  profiles  at  120  m  coinciding  with  the 
velocity  maximum.  The  kinetic  energy  of  the  central  Arctic  Ocean 
resides  almost  entirely  in  the  eddies  with  the  mean  currents  ac¬ 
counting  for  less  than  one  per  cent.  Although  detailed  surveys  of 
individual  eddies  are  lacking,  most  of  them  seem  to  have  a  diameter 
of  about  10  km.  On  that  basis,  it  is  estimated  that  eddies  cover 
between  10  and  20  per  cent  of  the  area  of  the  western  Arctic  Ocean. 
Their  sense  of  rotation  is  anticyclonic  or  clockwise.  Water  prop¬ 
erties  within  the  arctic  eddies  differ  from  those  without  and 
suggest  that  these  features  originate  in  the  shallow  waters  of  the 
Chukchi  Sea. 

A  number  of  possible  generating  mechanisms  for  these  eddies 
have  been  suggested  including  the  stress  of  wind  and  drifting  ice, 
brine  convection  accompanying  sea  ice  growth,  and  an  instability 
of  the  mean  baroclinic  current.  Of  these,  the  baroclinic  insta¬ 
bility  hypothesis  seems  most  plausible  and  is  the  only  one  examined 
here  in  detail. 

It  is  generally  accepted  that  the  migrating  cyclonic  systems 
of  the  atmosphere  as  well  as  the  turbulent  flow  in  certain  labora¬ 
tory  experiments  with  rotating  fluids  result  from  the  instability 
to  small  disturbances  of  the  mean  stratified  shear  flow.  This 
identification  of  baroclinic  instability  as  a  mechanism  for  gener¬ 
ating  eddies  in  the  atmosphere  and  laboratory  has  produced  a  body 
of  theory  which  has  been  summarized  and  reviewed  by  Charney  (1973) 
and  Kuo  (1973).  It  has  also  been  suggested  that  mesoscale  eddies 
observed  in  the  ocean  may  be  initiated  by  the  same  instability 


2 


mechanism  (Orlanski  and  Cox,  1973;  Gill  et  al,  1974). 

One  of  the  simplest  theoretical  demonstrations  of  baroclinic 
instability  is  that  by  Eady  (1949)  who  showed  that  for  a  basic 
state  in  which  the  vertical  shear,  VSis&lci  frequency  and  Coriolis 
parameter  are  all  constant,  certain  intermediate  wavelengths  are 
unstable  and  tend  to  grow.  Since  vertical  shear  and  VSisaia  fre¬ 
quency  are  not  actually  constant  in  the  ocean,  models  with  more 
realistic  structure  are  desirable  for  comparison  with  data. 

Eady's  theory  was  generalized  by  Williams  (1974)  to  include  mean 
vertical  profiles  which  vary  according  to  a  power  law.  He  further 
included  frictional  effects  by  introducing  Ekman  pumping  at  rigid 
upper  and  lower  boundaries  (Williams  and  Robinson,  1974) .  The 
results  were  successfully  applied  to  interpret  rotating  tank  ex¬ 
periments.  For  both  Eady's  and  Williams'  problems,  solutions  are 
in  terms  of  hyperbolic  functions. 

In  the  present  paper  Eady's  problem  is  extended  to  include 
mean  shear  and  VaisalS  frequency  which  vary  exponentially  with 
depth.  Analytical  solutions  to  the  problem  are  in  the  form  of  mod¬ 
ified  Bessel  functions.  Friction  is  included  in  the  form  of  Ekman 
pumping  at  either  rigid  or  free  surfaces.  A  model  similar  to  that 
developed,  in  this  paper  has  been  described  by  Hart  and  Killworth 
(1976)  who  also  applied  their  results  to  the  Arctic  Ocean.  The 
model  used  here  differs  by  including  free-surface  Ekman  pumping  and 
a  discussion  of  vertical  wave  structure. 

Observed  eddies  in  the  Atlantic  Ocean  are  in  the  neighborhood 
of  200  km  in  diameter  and  apparently  have  their  amplitude  maximum  at 
the  surface.  The  arctic  eddies  are  only  about  10  km  in  diameter  and 
have  their  maximum  amplitude  in  the  pycnocline  well  below  the 
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surface . 


These  differences  between  eddies  in  the  Atlantic  and 


Arctic  Oceans  suggested  an  investigation  of  factors  influencing  the 
horizontal  scale  and  vertical  profile  of  the  disturbance.  The 
shape  of  the  mean  shear  and  stability  profiles  as  well  as  the  pres¬ 
ence  of  the  ice  cover  on  the  Arctic  Ocean  are  expected  to  be  impor¬ 
tant  factors. 

2 .  The  Exponential  Eady  Problem. 

The  analysis  is  based  on  the  quasi-geostrophic  potential  vor- 
ticity  equation  in  a  form  used  extensively  in  oceanography  and 
meteorology.  For  an  ocean  with  a  mean  state  varying  only  with 
depth  and  which  is  disturbed  by  small  perturbations  in  pressure  and 
density  we  have 


and 


4- 


where  the  perturbations,  p*  and  f  ,  are  much  smaller  than  the  mean 
values.  If  there  is  a  geostrophic  balance  so  that 


u 


and 


/I  f  zy. 


*  x 


where 


is  the  geostrophic  stream  function 
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then  the  quasi-geostrophic  equation  may  be  written 


D 

D-b 


.  -h  -A-  / _f — 

v  T  ±7L\rJ- ^ 


o 

(2.1) 


where  directions  are  x,  east;  y,  north;  z,  vertically  upward;  and 
the  north  and  east  velocity  components  are  u  and  v.  The  Coriolis 
parameter,  f,  has  its  usual  meaning.  Other  symbols  are  defined  as 
follows : 


Total  derivative. 


D-b 


la 


4-  A r 


Horizontal  Laplacian, 


vaisaia  frequency. 


Beta  parameter. 


^  -z- 


The  quantity  in  brackets  is  the  potential  vorticity,  composed 
of  the  relative  vorticity  and  the  stretching  effect  of  density 
gradients  on  vortex  lines.  The  final  term  of  the  equation  repre¬ 
sents  the  effect  of  meridional  motions  on  potential  vorticity.  For 
the  derivation  of  this  equation  and  discussion,  see  Pedlosky  (1979). 

Assuming  a  stream  function  composed  of  steady  and  perturbation 


parts 


aj  V(z.)  + 
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where  the  steady  eastward  flow  varies  only  with  depth,  the  linear¬ 
ized  form  of  (2.1)  becomes 


( 


•  vy  + 


±{£±f 


)+PZt 


(2.2) 

This  equation  may  be  separated  by  introducing  a  perturbation 
stream  function  in  the  form  of  a  zonally-propagating  wave, 

ikU-c-t) 

f  =  FC2)  sLv\  e. 


to  give  the  differential  equation  for  vertical  structure 


+ 


1  (£  as)'. 

J2/ 


(2.3) 

2 

If  both  N  and  dU/dz  have  the  same  functional  form  and  6=0, 
the  final  term  in  brackets  vanishes  and  the  singularity  at  U  =  c  is 
circumvented.  This  has  been  called  the  generalized  Eady  problem  by 
Williams  (1974)  who  used  a  power-law  dependence  for  N"  and  dU/dz. 
Here  we  introduce  an  exponential  depth  dependence  for  those  quanti¬ 
ties  . 


Let  the  basic  density  variation  take  an  exponential  form 

/?(■=)  -  /f  -f  A/°  (/  -  £**) 


where  ps  is  the  density  value  at  the  surface  and  Ap  is  the  differ¬ 
ence  in  density  between  the  surface  and  great  depths. 
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The  VaisSla  frequency  in  this  case  becomes, 


d  A  p  & _  A  /° 

^  +  A  P  Q -  t.**)  ~  'Tf 


to  a  good  degree  of  approximation  since,  for  the  ocean, 

/°s  =  JO3  ki-m3,  ^  A/3  =  1  kj  ■  yyi'3 


and  the  scale  depth  of  the  pycnocline, 

o<  =  1  o  “3 


The  value  of  the  VSisSIS.  frequency  at  the  surface  is 


so  that 


The  basic  current  takes  the  same  exponential  form, 


c<  ZL 


These  profiles  have  a  resemblance  to  those  in  many  regions  of  the 
ocean  including  the  Arctic.  It  can  be  shown  that  the  association  of 
exponential  stratification  and  exponential  shear  implies  isopycnal 
surfaces  which  slope  uniformly.  Note  that  for  these  mean  conditions 
the  Richardson  number  increases  exponentially  with  depth, 
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Introducing  these  exponential  mean  conditions  into  (2.3)  produces 


the  non-dimensional  differential  equation, 

CL  -  ii-  ^  r  f  -o 

<J  z1  Jz. 

(2.4) 


Circumflexes  indicate  parameters  which  have  been  scaled  by  the  depth 
of  the  pycnocline,  z  =  oz,  and  by  an  internal  Rossby  radius  of  de¬ 
formation  , 


and 


Introduction  of  a  change  of  variable, 


F  (21)  = 

A 

where  A  71  /£, 

-  2  ^  ^ 


transforms  (2.4)  into  a  standard  form  of  the  Bessel  equation, 


_L  iZ. 

if  ^ 


-  0  +  p  =  ° 

(2.5) 


which  has  a  solution  in  terms  of  modified  Bessel  functions  of  the 
first  order, 
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where  A  and  B  are  arbitrary  constants.  Returning  to  the  original 
dependent  variable,  the  solution  becomes 


FC-fO  =  Af  I,(f)  -t-ZfK.Cf) 


(2.6) 


Boundary  conditions  at  the  surface  and  bottom  are  imposed  by 
Ekman  pumping  as  introduced  by  Charney  and  Eliassen  (1949) .  Vor- 
ticity  of  the  interior  frictionless  flow  leads  to  horizontal  mass 
divergence  in  the  Ekman  boundary  layer  and  consequent  vertical  flow 
at  the  interface  between  the  interior  and  the  boundary  layer.  Thus 
frictional  effects  in  the  boundary  layers  are  communicated  to  the 
interior  as  vertical  velocities.  The  Ekman  pumping  boundary  condi¬ 
tions  imposed  on  the  vertical  velocity. 


are  either 


U A 


for  a  rigid  upper  surface, 
such  as  an  ice  cover 


2. 


(2.7) 


or 


CM 


u 


for  a  stress-free  upper 
surface  as  in  the  case  of 
an  open  ocean 


(2.8) 
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For  a  lower  rigid  surface  the  Ekman  pumping  condition  is 


The  geostrophic  vorticity  is  defined  as 


-  vy  - 


T» 


he  scale  depth  of  the  frictional  boundary  layer  is 


De.  =  (2ii 


where  is  a  kinematic  eddy  viscosity. 

Substitution  of  (2.6)  into  the  boundary  conditions  gives  two 


equations  in  terms  of  the  two  unknown  amplitudes,  A  and  B.  For  rigid 


upper  and  lower  boundaries,  the  equations  are 


f  ( i  -  c)  y  J0(f .)  -  T,  (f .)  -iyy  L(fo) 


-+ 


£  [c  - 1 )  *  -  Wf  0  U  ^  fr*)_ 


=  o 


(2.9) 


and 


A\  (i*  -  ^ 


-hi  Y  y  I,(K 


_  r / a  -a 

y  3  (t-e. 


(2.10) 


where 


t 


r  EV 


and 


a  -  H /i- 
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while  for  a  free  upper  surface  (2.9)  is  replaced  by 

/\  [(/-£)  Y  z.(f.)  -  l(f.)  -  ^ 


+.  3  f(  £  -  0  y  X.  (f .)  -  K,(  f.)  + 


(2.11) 


where  H  is  the  dimensionless  ocean  depth,  K.  p  ,  and  the  dimension¬ 
less  phase  velocity  is  defined  by  c  =  c/U  .  The  Ekman  pumping 
parameters  are 

Y  =  De  Ns 

for  a  rigid  surface 

and 

v  =  K  •<  ^  An, 

for  a  free  surface. 

For  other  than  trivial  solutions,  the  determinant  of  the  co¬ 
efficients,  A  and  B,  must  vanish.  This  provides  a  dispersion  equa¬ 
tion  relating  c  and  y.  For  real  values  of  y,  c  is  generally  complex 

a 

with  the  real  part,  representing  the  phase  speed  and  the  imagi¬ 

nary  part,  c^,  representing  the  growth  or  decay  of  the  wave. 

Positive  values  of  c^  indicate  growing  waves  or  instability.  In 
this  case  the  dispersion  equation  is  a  quadratic  with  two  roots  and 
the  solutions  consist  of  either  a  pair  of  neutral  waves,  both  with 
Cj_  =  0,  or  one  growing  and  one  decaying  wave.  The  roots,  calculated 
with  the  aid  of  a  digital  computer,  are  plotted  in  Fig.  4  for  H  =  5 
and  in  Fig.  5  for  H  =  10.  Three  types  of  boundary  conditions  were 
used:  frictionless,  free-rigid  and  rigid-rigid.  The  phase  speeds, 

cr,  are  plotted  for  both  the  growing  and  decaying  waves.  Note  that 
for  free-rigid  boundary  conditions  the  lower  left-hand  branch  of 
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the  phase  speed  curve  is  associated  with  growing  waves  while  for 
the  rigid-rigid  case  it  is  the  upper  branch  which  corresponds  to 
growing  waves.  The  growth  rate,  yci ,  is  plotted  only  for  the 
growing  waves.  The  examples  all  show  a  maximum  growth  rate  at 
some  intermediate  wave  number  and  a  critical  wave  number  beyond 
which  no  growth  takes  place.  These  results  are  qualitatively  sim¬ 
ilar  to  results  by  others  for  Eady-type  models.  When  the  wave  is 
unchanging  in  the  y-direction,  m  =  0  and  y  =  k  to  give  the  highest 
values  of  the  wave  number  for  the  growth  maximum  and  for  the 

critical  point  beyond  which  growth  ceases.  The  scaled  water  depth 

*** 

in  Fig.  4,  H  =  5,  corresponds  generally  with  mid-ocean  conditions 
where  the  scale  depth  of  the  pycnocline  is  about  1  km  and  the  total 
water  depth  about  5  km. 

The  highest  growth  rate  is  found  for  the  case  with  no  Ekman 

A  A 

pumping  where  yc^  =  0.134  at  a  wave  number  of  0.8  when  H  =  5. 

Beyond  the  critical  wave  number,  f  =  1.18,  no  growth  occurs  and 
only  neutral  waves  exist.  When  Ekman  pumping  is  introduced  the 
growth  rate  is  reduced  and  the  wave  numbers  of  the  maximum  and 
critical  points  are  shifted  to  lower  values.  For  values  of  Y  =  0.5 
and  W  =0,05,  chosen  to  correspond  roughly  to  oceanic  conditions,  the 
dimensionless  growth  rate  is  reduced  to  0.068  for  a  free  upper  sur¬ 
face  and  to  0.02  for  a  rigid  upper  surface,  with  rigid  bottom  sur¬ 
faces  in  both  cases.  For  a  larger  ratio  of  pycnocline  depth  to 
total  depth,  H  =  10,  growth  rates  as  well  as  phase  speeds  are  less 
(Fig.  5)  . 

The  phase  speeds  of  the  decaying  and  growing  waves  coincide  in 
the  frictionless  case,  while  the  phase  speeds  of  the  neutral  waves 
split  into  two  separate  branches.  The  presence  of  Ekman  layers 
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splits  the  phase  speeds  for  the  growing  and  decaying  waves  also  into 
two  separate  branches. 


Vertical  profiles  for  the  fastest  growing  wave  are  shown  in 
Figs.  6 ,  7  and  8  for  three  different  boundary  conditions.  The  rela¬ 
tive  perturbation  stream  function  may  be  represented  as 


where  the  depth  function,  F,  has  been  decomposed  into  real  and 
imaginary  parts.  The  relative  magnitude  and  phase  of  the  perturba¬ 
tion  velocity  are  given  by 


and 


The  perturbation  velocity  at  the  upper  surface  is  greater  than 
that  on  the  lower  surface  in  all  of  these  cases,  although  in  Eady's 
original  problem  the  velocity  profile  was  symmetrical  about  the 
middle  level.  The  solutions  are  monotonic  functions  of  depth  in 
both  problems  and  there  are  no  maxima  at  intermediate  depths  either 
with  or  without  Ekman  layers. 

The  negative  phase  which  is  plotted  in  the  figures  presents  a 
lateral  view  of  wave  tilt  in  an  east-west  direction.  The  backward 
tilt  of  growing  waves  as  the  surface  is  approached  is  characteristic 
of  unstable  baroclinic  waves. 

Linear  perturbation  theory  predicts  conditions  only  at  the  on¬ 
set  of  instability  and  does  not  describe  later  developments  which 
occur  as  the  unstable  waves  continue  to  amplify.  However,  it  is 
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expected  that  the  wave  which  initially  grows  fastest  will  continue 
to  dominate  the  spectrum  of  unstable  waves  and  that  it  will  maintain 
the  general  character  of  its  original  profile. 

The  theory  presented  here  does  not  include  the  beta-effect  as 
an  approximation  to  the  latitudinal  variation  of  the  Coriolis 
parameter.  This  neglect  seems  justifiable  in  the  present  study 
since  it  has  been  shown  by  Green  (1960)  in  numerical  investigations 
that  although  the  inclusion  of  a  finite  beta-parameter  in  the  Eady 
problem  produces  an  additional  weaker  instability  maximum  at  low 
wave  numbers  and  also  destabilizes  the  high  wave  numbers,  it  does 
not  significantly  alter  either  the  wave  number  or  the  growth  rate 
of  the  fastest-growing  wave  which  remains  essentially  the  same  as 
in  Eady's  problem. 

3 .  Application  to  the  Arctic  Ocean. 

First  we  examine  the  possibility  that  the  arctic  eddies  are 
generated  in  offshore  locations  near  where  they  are  observed. 
Examples  of  profiles  for  geostrophic  current  and  vais&ia  frequency 
plotted  in  Figs.  9  and  10  show  that  an  exponential  curve  roughly 
fits  conditions  in  the  central  Arctic  Ocean.  Appropriate  constants 
based  on  these  data  are: 


The  Ekman  pumping  parameters  are  then 


and 


W  -  0.&5 ' 
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The  exponential  model  is  stable  for  Y  >  1.5  in  the  case  of  a 

<s 

rigid  top  representing  the  permanent  ice  cover  for  H  =  10.  Hence 
the  central  Arctic  Ocean  is  stable  to  small  disturbances  and  the 
eddies  cannot  originate  there  by  this  process.  If  we  choose  an 
improbably  low  value  for  the  Ekman  pumping  parameter  of  Y  =  0.5,  the 
flow  is  unstable  but  the  growth  rate 


corresponding  to  an  e-folding  time  of  612  days,  is  much  too  long  for 
significant  growth  in  this  situation.  The  half  wavelength  for  the 
fastest-growing  wave  in  this  case  would  be 


which  is  an  order  of  magnitude  larger  than  the  observed  diameter. 

If  the  arctic  eddies  are  formed  in  the  Arctic  Ocean  by  baro- 
clinic  instability,  they  must  originate  at  some  other  location 
where  conditions  are  less  stable.  At  the  edge  of  the  Alaskan  Con¬ 
tinental  shelf,  higher  mean  shear,  shallower  depth  and  lack  of  ice 
cover  during  late  summer  all  enhance  instability.  An  oceanographic 
section  made  by  the  USS  Staten  Island  in  September  1959  across  the 
Alaskan  shelf  north  of  Point  Barrow  (U.  S.  Naval  Oceanographic 
Office,  1963)  provided  guidance  in  selecting  parameters  for  this 
region:  Ng  =  10”^s”^,  LTS  =  10”^m«s“^,  =  10m,  a  =  10  “m 

FI  =  1,000m  (Y  =  0.5  and  W  =  0.05).  For  free  upper  and  rigid  lower 
boundries,  the  flow  is  unstable  with  a  scaled  maximum  growth  rate 
Of  0.055  at  a  scaled  wave  number  of  0.6  corresponding  to  a  dimen- 
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sional  growth  rate  of  15  days  with  a  half  wavelenght  of  37  km. 

This  rapid  growth  rate  indicates  that  mesoscale  eddies  should  form 
in  this  area.  The  dimensions  of  the  predicted  growing  disturbance 
are  roughly  comparable  to  observed  eddy  diameters.  Eddies  thus 
formed  would  be  advected  by  large-scale  currents  to  the  sites 
where  they  were  observed  in  the  deep  ocean  away  from  shelf  areas. 

On  the  other  side  of  the  Arctic  Ocean  in  the  area  northwest  of 
Svalbard  there  is  another  region  of  strong  shear  which  can  also  be 
examined  for  instability.  Warm  water  from  the  Atlantic  flows  north¬ 
ward  into  the  Arctic  Ocean  on  the  eastern  side  of  the  strait  be¬ 
tween  Greenland  and  Svalbard  while  cold  arctic  surface  water  leaves 
on  the  western  side.  Values  for  the  critical  parameters  in  that 
area  may  be  assigned  on  the  basis  of  a  survey  by  Newton  (1977): 

N  =  10”^s“l,  U  =  0.065m»s  ^ ,  D  =  7m,  a  =  1. 25xl0*"“m“^  and 

S  ^  0 

H  =  800m  so  that  y  =  0.5  and  w  =  0.05.  The  flow  is  unstable  with  an 
e-folding  time  of  18  days  for  the  fastest  growing  half  wavelength 
of  30  km.  These  results  are  similar  to  those  in  the  Alaskan  shelf 
region  but  it  is  not  likely  that  the  observed  eddies  originate  so 
far  away.  Oceanographic  surveys  are  being  planned  in  the  eastern 
Arctic  Ocean  which  should  provide  information  on  eddy  population 
north  of  Greenland  and  Svalbard.  Eddies  there  would  be  involved  in 
the  exchange  of  heat  and  salt  between  the  Atlantic  and  Arctic  Oceans 
and  may  be  a  factor  in  climatic'  change. 

4 .  Vorticity  diffusion  and  the  subsurface  current. 

Current  speeds  in  the  disturbances  generated  by  baroclinic  in¬ 
stability  have  their  maximum  at  the  top  pf  the  model  which  corres¬ 
ponds  to  the  base  of  the  Ekman  layer  at  a  depth  of  about  10  m.  Some 
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explanation  is  required  for  the  fact  that  the  observed  velocity 
maximum  is  much  deeper,  averaging  120  m.  If  the  eddies  are  indeed 
generated  along  the  edge  of  the  Alaskan  shelf  during  summer  when 
open  water  prevails  there,  then  the  current  maximum  will  be  near 
the  surface  at  the  base  of  the  Ekman  layer.  If  these  eddies  are 
then  carried  by  large-scale  currents  into  the  central  ice-covered 
ocean,  the  high  surface  velocities  will  decrease  rapidly  through 
friction  against  the  ice  leaving  a  subsurface  maximum.  This 
hypothesis  can  be  considered  quantitatively  as  a  problem  of  upward 
vorticity  diffusion  from  an  eddy  into  a  surface  layer  maintained 
at  zero  vorticity. 

First  we  need  to  examine  the  relation  between  the  Ekman  layer 
and  the  mixed  layer  which  extends  downward  to  a  depth  of  about 
50  m.  The  mixed  layer,  which  is  maintained  throughout  most  of  the 
year  by  brine  convection  resulting  from  ice  formation,  will  be 
spun  down  rapidly  through  friction  against  the  overlying  ice  cover. 
Secondary  circulations  will  cause  an  efficient  Ekman  spindown  of 
the  entire  layer  on  a  time  scale  of 

X  =  /■cxf” 

where  Hm  is  the  mixed  layer  depth  (Holton,  1972) .  Choosing  repre¬ 
sentative  values  for  the  Arctic  Ocean  of  Hm  =  40  m,  f  =  1.4  X  10“^s~^ 
and  'p*  =  7  X  10~^m2s“^,  we  find  a  rapid  response  time  of  16  hr. 

Below  the  mixed  layer  the  strong  stratification  will  suppress 
secondary  circulation  so  that  momentum  and  vorticity  will  be  lost 
upward  only  by  vertical  eddy  diffusion.  Then  the  time  scale  will  be 
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I 


where  H^is  taken  as  the  distance  from  the  base  of  the  mixed  layer  to 

the  subsurface  velocity  peak,  about  ICO  m,  and  must  have  a  re- 

-4  2 

duced  value  of,  say,  10  m“/s  so  that 


Thus  the  response  time  of  the  mixed  layer  at  the  surface  is  negli¬ 
gibly  small  in  comparison  with  that  of  the  stratified  layers  below 
it.  We  may  then  assume  that  the  frictional  effect  of  the  ice  is 
transmitted  directly  to  the  base  of  the  mixed  layer. 

We  examine  the  case  of  an  initial  vorticity  profile  with  a 
maximum  at  the  top,  that  is  at  the  base  of  the  mixed  layer,  where 
the  vorticity  remains  zero  due  to  friction  and  trace  the  changes  in 
the  profile  with  time  as  a  subsurface  maximum  develops  which  resem¬ 
bles  the  observations.  Since  the  horizontal  and  vertical  structures 
are  separable  and  since  vertical  vorticity  is  a  horizontal  deriva¬ 
tive  of  velocity,  the  vertical  profile  derived  for  vorticity  will 
also  apply  to  velocity. 

In  the  absence  of  vertical  stretching  and  horizontal  divergence, 
the  linearized  equation  for  vertical  vorticity  is  a  diffusion  equa¬ 
tion. 


o 


(4.1) 


where  vorticity 


^  at  _  ^  u. 
d  x  } 
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In  terms  of  a  Fourier  integral,  (4.1)  has  a  solution 

/  _  (z-z'/Aja  t 

z.ty  =  — _ 


J  z' 


(4.2) 


where  the  initial  vorticity  is 


f (z‘)  . 

4  ^ 


and 


a  S 


The  boundary  conditions  are 

Z-  -  O  for-  t-  0 

Z_  —  oO  elH  ^  _ 


Although  the  vertical  profile  derived  from  the  previous  baroclinic 
instability  problem  has  the  form  of  a  modified  Bessel  function,  we 
approximate  that  shape  here  with  an  exponential.  Assume  an  expo¬ 
nential  profile  which  has  a  mirror  image  above  to  maintain  zero 
vorticity  at  z  =  0.  , 


-f/  k-  2:  >  o 
■(  .  r  <  O 


(x-x')W^fc 


oO 

c 
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o 


or 


z.t)  -  - 


A 


J 


Z(i?»t.) 


‘A 


ft2 

c  /c 


CK 


-cO 


eO 


+  t 


a-  /  -Z?'  -  (  Z  -  3  Z  £*)]*//-&-£ 


Jz' 


(4.3) 


In  the  first  integral  of  (4.3),  let 


'  -  (z  tcQ 


2  (>  t; 5/1 


=  e 


and  in  the  second  integral,  let 

Z  '  -  (Z  -  2  a) 

2  0i-VA 


This  gives 

J  (X*)  =  /)  77-"^- 


2  -f  2Ut«V 

e:6Ae 
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Introducing  the  error  function  notation, 


e-rf 


we  have 


V. 


4  t 


2£U-2^\ 

2  o>  t)'4 y 


(4.4) 


Equation  (4.4) 
for  vorticity, 


may  be  evaluated  in 
;  depth,  az; 


terms  of  normalized 

// 

and  time,  C\  *  ^ 


variables  : 
(Fig.  11) . 


Vorticity  diffuses  upward  and  downward  so  that  a  subsurface  maximum 
develops  which  both  decreases  with  time  and  descends  to  deeper 
levels.  These  theoretical  profiles  closely  resemble  the  observed 
profiles,  for  example  see  Fig.  1. 

From  the  previous  results  on  baroclinic  instability  we  may 
assume  that  1/a  =  200  m.  The  choice  of  an  eddy  coefficient  is  more 
difficult  but  experience  in  steep  density  gradients  of  other  oceans 
suggests  a  value  of  ^  =  lO'^m^/s  for  an  order  of  magnitude  cal¬ 

culation.  In  the  Arctic  Ocean  the  velocity  maximum  occurs  at  100 
to  150  m.  The  mixed  layer  depth  of  50  m  must  be  subtracted  to  cor¬ 
respond  to  the  theoretical  problem.  Then  az  =  0.25  to  0.5. 
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Referring  to  Fig.  11  it  is  seen  that  this  corresponds  to  a  dimen¬ 
sional  decay  time  ranging  from  46  to  135  days,  indicating  that  the 
observed  eddies  had  spent  from  1  to  6  months  beneath  the  ice.  The 
decay  time  is  directly  proportional  to  choice  of  the  eddy  coeffi¬ 
cient  value  which  is  admittedly  a  poorly  known  quantity  and  this 
eddy  age  can  only  be  considered  suggestive. 

5 .  Summary  and  discussion. 

The  vigorous  mesoscale  features  which  are  so  abundant  in  the 
Arctic  Ocean  require  a  source  of  energy  for  their  generation  and 
maintenance.  Although  wind  energy  has  been  suggested  as  a  source, 
the  lack  of  correlation  between  winds  and  mesoscale  motion  makes 
this  an  unlikely  relationship.  Furthermore  the  arctic  eddies  are 
much  smaller  in  horizontal  scale  than  synoptic  weather  systems 
although  a  similar  size  might  be  expected  if  the  eddies  were 
driven  by  winds.  Another  possible  energy  source  which  has  been 
suggested  is  thermodynamic  in  nature  and  associated  with  the  freez¬ 
ing  process.  During  winter,  areas  of  open  water  appear  in  the  ice 
pack  through  either  breakup  or  differential  motion  of  ice  floes. 

Ice  grows  quickly  in  these  open  areas  at  the  extremely  low  air 
temperatures  which  prevail.  As  the  ice  grows,  salt  is  expelled  and 
heavy  brine  sinks  through  the  mixed  layer  causing  intense  local 
convection.  But  it  is  difficult  to  understand  how  this  shallow 
convection  can  produce  the  observed  redistribution  of  density  at 
the  much  deeper  levels  where  eddies  are  found.  Also  there  is  here 
again  a  mismatch  in  scales  between  the  sizes  of  source  and  eddy. 
Open  areas  where  rapid  freezing  occurs  are  usually  tens  to  hundreds 
of  meters  across  while  the  eddies  are  about  10  kilometers  in 
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diameter.  A  more  plausible  source  for  eddy  kinetic  energy  than 
either  of  these  suggestions  is  the  available  potential  energy  re¬ 
siding  in  the  inclined  density  surfaces  accompanying  the  mean 
geostrophic  circulation  of  the  Arctic  Ocean.  The  mean  vertical 
shear  has  been  shown  to  be  unstable  in  the  Barrow  Current  north  of 
Alaska  but  not  in  the  central  Arctic  Ocean  where  eddies  were  ob¬ 
served.  The  small-amplitude  theory  employed  predicts  an  eddy  size 
close  to  that  observed  and  it  predicts  a  growth  rate  fast  enough 
to  indicate  rapid  formation  of  these  features.  Although  the  theory 
cannot  describe  the  evolution  of  an  eddy,  it  is  assumed  that  the 
unstable  waves  continue  to  develop  until  meanders  form  which 
eventually  cut  off  to  become  isolated  eddies. 

The  baroclinic  instability  theory  used  here  applies  to  a  ver¬ 
tical  shear  only  and  is  appropriate  to  mid-ocean  circulation  where 
the  horizontal  shear  may  be  considered  negligible.  It  incorporates 
an  exponential  mean  shear  and  stratification  as  well  as  Ekman 
pumping  in  an  attempt  to  more  closely  fit  actual  oceanic  conditions. 
The  conclusion  that  the  central  Arctic  Ocean  is  stable  and  that 
eddies  do  not  originate  there  by  this  process  seems  quite  definite. 
They  must  be  formed  elsewhere,  apparently  in  the  Barrow  Current 
and  application  of  the  theory  to  this  current  confirms  that  it 
could  be  the  source. 

There  are  two  assumptions  inherent  in  the  theory  which  are  not 
exactly  met  in  the  Barrow  Current  however.  The  Current  is  a  rela¬ 
tively  narrow  jet  which  has  significant  horizontal  shear  and  it 
occurs  along  the  shelf  edge  where  the  topography  is  not  flat.  But 
it  is  not  believed  that  the  conclusions  reached  would  be  seriously 
altered  by  inclusion  of  a  horizontal  mean  shear  effects  and  depth 
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variations.  Although  no  quantitative  modeling  demonstrates  this, 
some  qualitative  arguments  make  the  omission  of  these  effects  ap¬ 
pear  relatively  unimportant. 

A  mean  horizontal  shear  may  lead  to  instability  and  it  is 
quite  possible  that  the  Barrow  Current  is  barotropically  as  well 
as  baroclinically  unstable.  The  combined  barotropic  and  baroclinic 
instability  problem  is  an  inseparable  one  and  solutions  to  it  have 
been  primarily  sought  by  numerical  methods.  Barotropic  instability 
could  only  be  expected  to  enhance  eddy  production.  Its  effect  on 
growth  rate  and  eddy  size  in  this  case  will  require  further  model¬ 
ing  in  the  future. 

Bottom  topography  such  as  the  continental  margin  in  the 
vicinity  of  the  Barrow  Current  is  a  departure  from  the  level  bottom 
assumed.  Theoretical  studies  of  topographic  effects  on  the  baro¬ 
clinic  instability  problem  cited  in  McWilliams  (1979)  review  show 
that  the  classical  Eady  waves  studied  here  are  not  strongly  influ¬ 
enced  by  depth  variations  and  that  instability  may  even  be  enhanced 
by  such  variations.  It  is  reasonable  that  bottom  topography  would 
not  seriously  modify  the  waves  studied  here  wThich  are  trapped  pri¬ 
marily  against  the  upper  surface  of  the  ocean. 
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ILLUSTRATIONS 


Fig.  1  -  Profile  of  horizontal  current  through  an  Arctic  Ocean 

eddy.  Solid  line  is  current  speed  and  dashed  line  is  direc¬ 
tion,  5  May  1975,  76°23'N  149°48'W  (Manley,  1981). 

Fig.  2  -  Profile  of  temperature  (T)  ,  salinity  (S)  ,  and  density  (<37^) 
for  the  same  eddy  shown  in  Fig.  1.  Solid  lines  are  for  condi¬ 
tions  outside  eddy  in  surrounding  waters  and  dashed  lines  are 
for  conditions  inside  eddy  (Manley,  1981) . 

Fig.  3  -  Locations  of  arctic  eddies  observed  during  drift  of  four 
AIDJEX  ice  camps.  Of  the  98  eddies  for  which  rotation  sense 
was  determinable,  95  were  anticyclonic  or  clockwise.  The 
three  cyclonic  eddies  are  circled. 

Fig.  4  -  Phase  speed  and  growth  rate  of  baroclinic  waves.  H  =  5. 
Solid  line  is  for  rigid  upper  and  lower  surfaces  with  no 
friction.  Dashed  line  is  for  Ekman  pumping  at  a  free  upper 
surface  and  rigid  lower  surface,  W  =0.05,  Y  =0.5.  Dotted 
line  is  for  Ekman  pumping  at  upper  and  lower  rigid  surfaces, 

Y  =  1. 

A 

Fig.  5  -  Phase  speed  and  growth  rate  of  baroclinic  waves.  H  =  10. 
Symbols  as  in  Fig.  4. 

Fig.  6  -  Relative  amplitude,  Mf  and  negative  phase, — ^  ,  of  fast¬ 

est-growing  wave.  Rigid  upper  and  lower  surfaces  with  no 
Ekman  pumping.  H  =  5  and  R  =  0.8. 

Fig.  7  -  Relative  amplitude,  M,  and  negative  phase,  ^  ,  of  fast¬ 

est-growing  wave.  Ekman  pumping  with  free  upper  and  rigid 
lower  surface.  ft  =  5  and  ft  =  0.7. 

Fig.  8  -  Relative  amplitude,  M,  and  negative  phase,  ~~  Q  ,  of  fast¬ 
est-growing  wave.  Ekman  pumping  with  rigid  upper  and  lower 
surfaces.  H  =  5  and  ft  =  0.4. 

Fig.  9  -  Geostrophic  currents  in  the  Arctic  Ocean.  Profiles  (1) 

and  (2)  are  in  the  central  Arctic  Ocean,  (3)  is  near  the  shelf 
edge  north  of  Alaska  (Kusunoki,  1962). 

Fig.  10  -  Profiles  of  VSisSlS  frequency  in  the  central  Arctic  Ocean. 

Fig.  11  -  Development  of  a  vorticity  profile  with  a  subsurface  max¬ 
imum  from  an  initial  exponential  profile.  Vorticity  at  the 
upper  surface  is  kept  at  zero  by  frictional  dissipation. 


28 


Fig 


1 


TRUE  DIRECTION 

O  90  1B0  270  360 


29 


(SIGMfl-T) 


Fig.  2 


0- 

Q_ 


-h 

Q 

a 


+ 


o 

O 

o 

O 

CM 

cn 

4- 

o 

o 


-4- 

o 

o 

in 


-h 

o 

o 

CO 


■+■ 

a 

o 

r~ 


4 


30 


800 


78-  nlfj*''  M 9W  1 45W 


Fig.  3 


31 


Growth  Rate,  o  C:  Phase  Speed, 


Fig. 


4 


32 


0.15 


<cT 

<^o 

%  0.10 

a 

45 

4J 

o  0.05 

a 


0 


Fig.  6 


M  0  0.5  1.0  1.5 


34 


Fig.  7 


M 

-a 


35 


A 

Z 


5 


CURRENT  VELOCITY,  cm  /sec 

5  10  15 


37 


0 

50 

100 

150 

200 

250 

300 

350 


Fig.  10 


* 


.01 

- r~ 


N,  rad /sec 

.02  .03 


.04 

“I 


38 


DEPTH, az 


r  i  '  . 


11 


RELATIVE  VORTICIT Y ,  £/A 


39 


MANDATORY 
DISTRIBUTION  LIST 


FOR  UNCLASSIFIED  TECHNICAL  REPORTS,  REPRINTS,  &  FINAL  REPORTS 

PUBLISHED  BY  OCEANOGRAPHIC  CONTRACTORS 
OF  THE  OCEAN  SCIENCE  AND  TECHNOLOGY  DIVISION 
OF  THE  OFFICE  OF  NAVAL  RESEARCH 


1  Director  of  Defense  Research  and 
Engineering 

Office  of  the  Secretary  of  Defense 

Washington,  D.  C.  20301 

ATTN:  Office,  Assistant  Director 

( Research) 


1 

1 

1 

3 


Office  of  Naval  Research 
Arlington,  VA  22217 
ATTN:  (Code  102-C) 

ATTN:  (Code  200) 

ATTN:  (Code  460) 

ATTN:  (Code  480) 


6  Director 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 
ATTN:  Library,  Code  2620 


1  U.  S.  Naval  Research  Laboratory 
Code  2627 

Washington,  D.  C.  20375 


2  Office  of  Naval  Research  -  N.Y. 
715  Broadway 
New  York,  N.  Y.  10003 


12  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 


1  Commander 

Naval  Oceanographic  Office 
NSTL  Station 
Bay  St.  Louis,  MS  39522 
ATTN:  Code  02 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

'•  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

CU-2-81 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

ARCTIC  OCEAN  EDDIES  AND 
BAROCLINIC  INSTABILITY 

5.  TYPE  OF  REPORT  6  PERIOD  COVERED 

Technical  Report 

March  1981  to  May  1981 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author/*; 

Kenneth  Hunkins 

a.  CONTRACT  OR  GRANT  NUMBER/*; 

N00014-76-C-0004 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Lamont— Doherty  Geological  Observatory  of 
Columbia  University 

Palisades,  N.  Y.  10964 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

NR307-359 

1  1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

1  Dept,  of  Navy,  Office  of  Naval  Research 
Code  461,  Arlington,  VA.  22217 

12.  REPORT  DATE 

July  1981 

13.  NUMBER  OF  PAGES 

39 

14.  MONITORING  AGENCY  NAME  &  ADDRESS//*  different  from  Controlling  Office) 

Same  As  Above 

15.  SECURITY  CLASS,  (of  this  report) 

Unclassified 

1 5a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose 
of  the  United  States  Government. 

17.  DISTRIBUTION  STATEMENT  (of  the  abetract  entered  in  Block  20,  if  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  reverse  aide  if  neceesary  and  identify  by  block  number) 

Eddies,  Arctic  Ocean,  baroclinic  instability, 
mesoscale  variability 

20.  ABSTRACT  (Continue  on  reverae  aide  it  neceaaary  and  identify  by  block  number) 

Baroclinic  eddies  with  diameters  of  about  10  km  and  maximum 
current  speeds  of  about  50  cm/s  have  beeb  widely  observed  in  the 
central  Arctic  Ocean  north  of  Alaska  and  Canada.  The  possible 
origin  of  these  eddies  through  an  instability  of  th  mean  baro¬ 
clinic  flow  is  investigated  using  an  model  with  exponential 
profiles  of  mean  shear  and  Vaisala  frequency.  The  model  includes 
Ekman  pumping  at  a  rigid  bottom  and  at  either  a  free  or  rigid 

1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 

S/N  01 02-LF-01 4-6601 


DD  FORM 

1  J  AN  73 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (Whan  Data  Entered) 


upper  surface.  The  central  Arctic  Ocean  where  the  eddies  were 
found  is  baroclinically  stable  with  no  possibility  of  eddy  pro¬ 
duction.  If  eddies  are  spawned  by  this  mechanism,  they  must 
be  formed  at  a  site  far  from  where  they  are  observed.  On  the 
periphery  of  the  Arctic  Ocean  north  of  Alaska  the  combination  of 
greater  current  shear,  shallower  depth  and  lack  of  ice  cover 
leads  to  unstable  conditions  and  the  eddies  apparently  originate 
in  that  region.  The  instability  theory  predicts  maximum  velocity 
at  the  surface  instead  of  below  the  surface  as  observed*  Appar¬ 
ently  after  formation  in  the  open  water  the  eddies  are  advected 
beneath  the  ice  cover  and  dissipate  the  momentum  of  their  upper 
layers  against  the  ice.  This  is  demonstrated  by  calculations  for 
the  diffusion  of  vorticity  against  the  ice  in  the  case  of  an 
initial  exponential  profile.  A  subsurface  maximum  then  develops 
which  resembles  the  observed  profiles. 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWh*n  Data  Entered) 


COLUMBIA  LIBRARIES  OFFSITE 


90646606 


\ 


